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The collision of two gravitationally interacting, ultra-relativistic, extended sources is being examined. 
This investigation classifies the transverse distributions that are collided for fixed collision energy, 
according to whether one or two (a small and a large) apparent horizons may or may not be formed 
in a flat background in 4 dimensions. The study extends to the thermodynamical properties of 
the objects that are created, which exhibit a universal behavior in their entropy, and, suggests 
the elimination of the possibility in observing black holes (BHs) at the LHC in the absence of 
extra dimensions. On the other hand, including extra dimensions, and assuming that the matter is 
localized (dense) enough in those directions, opens new avenues in creating BHs at energies of the 
order of TeV. The investigation is carried further to AdSc, backgrounds and makes connections with 
the implications for the quark-gluon plasma (QGP) formation in heavy ion collisions. In particular, 
classes of the geometries found suggest that a BH is formed if and only if the (central collision) 
energy is sufficiently large compared to the transverse scale of the corresponding gauge theory side 
stress-tensor. This implies that when the scattering in the gravity description is mapped onto a 
heavy ion collision problem yields a result, which is in accordance with the current intuition and 
data: QGP is formed only at high enough energies compared to Kqcd, even for central processes. 
Incorporating weak coupling physics and in particular the Color Glass Condensate (CGC) model, a 
satisfactory fitting with the RHIC and the LHC data for multiplicities may be established. 



I. INTRODUCTION 

Gravitational ultra-relativistic collisions of two shock- 
waves generated by point-like sources on the transverse 
plane have been discussed in several contexts, [l|-[l8|. 
Earliest studies have shown that such collisions typ- 
ically create a single black hole (BH) given sufficiently 
large energy compared to the impact parameter involved. 
In particular, for central collisions, these and some fol- 
low up works [3, [l9l - [23| , have shown that a BH is formed 
for any (small) collision energy. Works [13, H^l general- 
ized the analysis to extended distributions and found the 
conditions under which a single BH may be formed. 

In this work further steps are taken in that direction. 
We classify the energy-momentum (transverse) distribu- 
tions that are collided according to (i) always a trapped 
horizon^ is formed independently on the amount of the 
collision energy (as in 14]), (ii) one trapped horizon is 
formed for sufficiently large energy (as in 10]) (iii) two 
(and sometimes more) trapped horizons, one small and 
one large, are formed for sufficiently large energy. The 
classification is done in the two cases of flat backgrounds 
and AdSs backgrounds where the creation of mini BHs 
(see also ^^) in the accelerators and the implications 
in heavy ion collisions are discussed respectively. We 
restrict the classification to transversally symmetric ge- 
ometries ^ in the spirit of trapped surfaces. Our analysis 
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^ Unless otherwise specified, we will use the terminology trapped 
horizon, trapped surface and BH interchangeably. 

^ The transverse symmetry will be defined rigorously for each back- 
ground. Other than this symmetry and three mild and physical 
assumptions, the classification is completely generic. 



assumes a zero impact parameter limit as this case suf- 
fices in order to illustrate our claims. 

BHs from shock-wave collisions is not a new concept. 
What is new here is: 

1. Category (iii) of trapped surfaces formed under a 
dynamical process has not been studied so far in the lit- 
erature and is one of the main results of our work. In 
fact, evidence of creation of small (in addition to large) 
BHs in flat backgrounds has not been reported up to to- 
day. The underlying set-ups consist of examples of the 
(known) non-uniqueness of the solutions regarding the 
trapped surface equations (see [l^). 

2. The small BHs seem to have a curious property: 
their size (and hence their entropy) decreases as the en- 
ergy increases and this behavior is exactly opposite to 
that of the large BHs. If this was an equilibrium state^ 
it would imply a negative temperature (T). There are 
a number of interesting scenarios to be discussed. An 
attractive explanation is that the resulting trapped sur- 
face, having two (apparent) horizons, will evolve into a 
Reissner-Nordstrom-like (RN) of BH whose inner and 
outer horizon share a similar behavior, as the energy 
varies, with the corresponding small and large apparent 
horizons that we flnd in our analysis. In fact, there ex- 
ists a critical energy such that the two surfaces have the 
same size, which by analogy to the RN BH, it could be 
the extremal limit. Lowering the energy further, in our 
case we find that no trapped horizons may be formed. In 
the RN case, it yields a naked singularity, which violates 



In fact it is not; this is the beginning of a very violent process. 
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the "cosmic censorship" conjecture. The role of the cor- 
responding scale of the electric field in the RN solution, is 
(possibly) played by the transverse width of the colliding 
distributions (see also [2^[13|)- In particular, working in 
an AdSs background, such a configuration could model 
[l9l - l23l . [28I - I30I I heavy ion collisions. As the inner horizon 
has a different thermodynamical behavior than the outer 
region, an RN-like scenario could possibly describe a hot- 
ter fireball in the center of the produced plasma, which 
absorbs heat from the surrounding shell, and a thermo- 
dynamically unstable periphery [3ll-[34 1 that undergoes a 
hadronic phase transition. 

3. In the high energy (E) limit, the entropy (S) of the 
(large) BHs for all the three classes of the colliding distri- 
butions and for both, flat and AdS backgrounds, exhibits 
a universal behavior. The entropy as a function of the 
cm. energy s = AE'^ and the transverse scale k of the 
distribution behaves as 5 ^ k^s and 5* ~ (s/fc^)^/'^ for 
flat and AdSs backgrounds respectively independently on 
the details of the profiles of the (bulk) distributions that 
are collided. In general, k can be s dependent (e.g. see 
(PO)) ). In fact, this universal behavior is independent on 
the impact parameter and any extra dimensions (if these 
ingredients are introduced) as long as the energy is suffi- 
ciently large. Applying the AdS/CFT results of this work 
in heavy ion collisions and incorporating weak coupling 
physics, a satisfactory fitting with the total multiplicities 
data may be achieved. This can be seen in fig. |S1 

4. Under the assumption that extra dimensions do not 
exist, this work motivates that the LHC is a "black hole 
safe" machine. The parameter that controls the BH for- 
mation is the dimensionless combination Ek/M^ > w 
where E is the collision energy, k the is transverse scale 
associated with the colliding distribution p and can be 
taken of the order of Aqcd, Mp is the Planck mass and 
w is a dimensionless number that is predicted from the 
trapped surface equations and naturally depends on p. 
Usually w is of the order of 1 or less. In particular, 
for the boosted Woods-Saxon profile for nuclear distri- 
butions, which is a suitable phenomenological model for 
collisions at LHC, w « 0.4. This implies that the colli- 
sion energy that is required to create BHs at the LHC 
is way beyond the machine's possibilities. Including ex- 
tra dimensions on the other hand, could change things 
dramatically. Towards the end of the conclusions, it is 
argued that BHs could be created even at TeV scales 
provided that the distributions that are collided are suf- 
ficiently dense in the compact directions. In other words, 
the underlying assumption is that the matter is confined 
almost completely on the 3-brane. 

5. Classes (ii) and (iii) for AdS backgrounds, advo- 
cate the idea that not any distribution of matter at low 
energy when collided, yields a BH. M apping this idea, 
through the AdS/CFTjsidictionary to heavy 
ion collisions as in [H-illai-liil implies that QGP ^ll 
is formed provided that the E of the colliding (dual pic- 



ture) glue-balls is sufficient. In particular, according to 
our calculations, E must be much larger than the trans- 
verse scale of the corresponding gauge-theory stress ten- 
sor (that can be taken of the order of Aqcd), even for 
a zero impact parameter collision. This is a totally new 
result and extends the results of the existing literature, 
which state that QGP, under a holographic investigation, 
is always formed for any small collision energy as long as 
the collision is (perfectly) central . In particular, as we 
show for AdS backgrounds, the parameter controlling the 
BH formation appears in the combination E/k > w (con- 
trast with point 4. above for flat backgrounds) where w 
depends on the particular (bulk, 5-dimensional) stress- 
tensor distribution. Central nuclei collisions at energies 
of a few eV's should not result to a QGP formation; this 
has been the initial motivation of this part of our inves- 
tigation. 

We organize this paper as follows. 

In section [H] we set up the problem and we briefly 
review the method that we will use in order to support 
our claims. 

Section IIIII deals with the classification of the energy 
momentum distributions for fiat backgrounds. Three 
classes of distributions are found and one example for 
each one is given. In addition, the thermodynamical con- 
sequences, the final product of the collision, as well as the 
implications to collisions at the LHC are discussed. 

Next section, section ITVl generalizes the discussion of 
section iHll to AdS backgrounds and gives particular em- 
phasis to the consequences in the applications of the 
gauge/gravity duality in heavy ion collisions. 

Finally, in section |V] we summarize our findings and 
suggest possible extensions of this analysis. In particu- 
lar, we include an extended discussion for extra fiat di- 
mensions and argue that the inclusion of such compact 
directions could make BHs creation at the LHC feasible. 



II. SETTING UP THE PROBLEM 

In this section we review very briefiy the trapped sur- 
face method that is required for our subsequent analysis. 
More details can be found in [H [ill2|j4i| while ex- 
tended notes exist in the appendices of [23|. 

We begin with the shock metric in a flat background 



created by a fast-moving source moving along the 



direction 



—2dx dx 
{x\x^) 



{x^,x^)5{x'^)dx 



+ 2 



■ dx 



(1) 



where (j) is the transverse profile of the shock and has 
dimensions of length, a;^ are light-cone coordinates (de- 
fined as x^ = [x^ ^x^) / ^/2) and i5 is the Dirac's function. 
This ansatz satisfies trivially all the Einstein's equations 
except the (-I--I-) component, which yields 
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^8ttGT++ = EGpd{x+), T++ = Ep5{x+) (2) 

allowing one to identify the stress energy tensor com- 
ponent r++ with the transverse Laplacian of </> (times 
5{x'^)/GY. In equation ^ we find it convenient to de- 
fine T++ by pulling out its total energy E and its delta 
function dependence; this defines the density p = p{x±). 

The trapped surface consists of two pieces, S+ and 
These are parametrized with the help of two functions, 
V'+ and V"- , which satisfy the following differential equa- 
tion 



vi(v± 



= 0. 



(3) 



Obviously, V^(/)± provides a source term for 'S/\ip±. The 
missing ingredient is the boundary conditions. Dropping 
the indices ± from tp± from now on assuming a zero 
impact parameter and identical shocks we have ■;/'+ = 
'0_ = ip. The boundary conditions then read 



c 



[v=^'V'V,.v] 



1=1.2 



(4) 



for some curve C, which defines the boundary of the 
trapped surface (so the curve C is the trapped horizon) 
and where both, S+ = S and S- = S end. The set of 
equations, equations ^ and dH) that define the trapped 
surface are generalized for AdS spaces in a straightfor- 
ward manner. In particular, the differential equation ([3]) 
for AdSs is replaced by 



where 0± is the scalar operator in AdS^ while the (left 
hand side of the) second boundary condition in (U) should 
be contracted with the metric tensor of AdSa. 



III. FLAT BACKGROUND BHS PRODUCTION 



= 8 => 



1. 



(7) 



Now taking into account that equation Vi0,, ^ EGp 
(see equation ^) implies l/x±{x_\_(j)')' ^ EGp, then 
equation ([7]) yields 



/p^^ xp{x, k)dx 



1 

'GE 



> 



f^^ xp{x, k)dx 



kx'^t 



GEk 



(8) 



where k a transverse scale whose inverse determines the 
width of the distribution p. Equation ([8]) is the basic 
equation of the flat background investigation and its con- 
sequence can be fit in three classes depending on p. The 
only assumptions for p is that it must be (i) everywhere 
positive (in space-time) in order to be a physical energy 
distribution^ , (ii) integrable in order to yield a finite total 
collision energy and (iii) it must have a monotonic fall-off 
in the sense d/dx{xp{x)) — has only a single real root*". 

In particular, assumption (ii) implies that the numer- 
ator in the left hand side of ([S]) is finite as x]^ — > oo. 



Assumptions (i) and (ii) combined, imply that the roots 
of the left hand side of ([5]) is at a;^ = oo and possibly 
at x'j_ — depending on the behavior of p at the origin. 
As it will be seen shortly, the behavior of p{x±) near 
x±_ = partitions the set of positive and integrable p's 
into three disconnected classes according to the details of 
the trapped horizons. 

The physical meaning of formula (|5]) is that given a 
total collision energy E, a trapped horizon exists if and 
only if there exists an x]_ such that the energy that is 
(5) enclosed inside the circle of radius a;]^ is equal (modulus 
constant pre-factors) to x'j^/G. An observation for (|S]) 
comes about when notting that the left hand side of the 
equation (after the => symbol) is a dimensionless quan- 
tity of the variable kxc- This implies that the dimen- 
sionless quantity that determines the existence and the 
behavior of the trapped surface(s) in a 4-dimensional fiat 
backgrounds is set by the product of two energy scales as 
E X k in units of G. This fact should be contrasted with 



General classification 



The solution to the system of equations and (|H) is 

i^{x^) = cPix^)-(f>{xl) (6) 

where the curve C : x± — x'j^ defines the trapped horizon, 
which is a 2-dimensional disk. The second boundary con- 
dition then specifies x']_ in terms of the data of the shock 
yielding^ 



* From now on we drop the dimensionless proportionality factors 
for simplicity as they do not alter qualitatively the final results. 
^ All the positive constants will be set equal to one from now on. 



This is the weak (and since Tj^ = it is also the strong) energy 
condition: T^'^UfiTii, = T ('^— )^ > for any time-like (and 

space-like) vector if and only if T++ = T ~ p > 0. 

By assumption as one expects that physical distributions either 
will increase, reach a maximum value and then decrease at in- 
finity or they will start from a maximum value (which can be 
infinite but integrable at the origin) and keep decreasing. Mul- 
tiplying by the measure x (of polar coordinates as the system is 
effectively 2-dimensional), turns either type of the distribution 
to a bump function, which vanishes at zero and at infinity and 
it has a maximum between. If this assumption is removed, then 
as we will see, more than two trapped horizons are allowed as 
the number of trapped surfaces increases with the number of the 
extrema of the left hand side of JSjl. In particular, assumption 
(iii) restricts the number of extrema to one or less (this is not a 
straightforward fact-see footnote in subsection IIII A"3l l . 
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FIG. 1: The function (r(l/4) - r(l/4, y^)/?/ as a function of 
y = kx±. The dashed line corresponds to a given energy E. 
This is an example of a distribution that always results to 
a single trapped horizon regardless from the (small) value of 
GkE. 

the next section, which deals with AdSs backgrounds. 
There, as will see, the behavior of the corresponding di- 
mensionless coupling is dramatically different and inter- 
esting. 
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FIG. 2: The function Erf(j/)/j/ as a function of j/ = kx±. The 
dashed lines correspond to lines of constant energy E. This is 
an example of a distribution that results to a single trapped 
horizon but only for sufficiently small values of 1/GkE and 
hence sufficiently large values of E in units of 1/kG. In par- 
ticular, the top dashed line does not correspond to a trapped 
horizon formation while the bottom line does. 



2. Marginal case: one black hole for sufficiently large 
energy 



1. Always a single black hole for any energy 

This case corresponds to the distributions, which 
near the origin behave as p{x±_) ^ -I- (sub- leading) 
where 2 > n > 1 ^ ^. This implies that the left hand 
side of ([5]) diverges at = and goes to zero as 
x']_ — )■ oo and as a result, there is always a single^*^ 
solution for any E to equation ([8]) for some x']^. As 
an example, we consider the (everywhere positive and 
integrable) distribution p ~ k-^e~'' /{kxj_)^/^ , 
which (for completeness) results to the shock 
- £;Gyfcil2F2({l/4,l/4};{5/4,5/4};-(fcxi)2) 
where pFq the generalized hypergeometric function. 
Here k is some transverse scale that fixes the width 
of the distribution. The left hand side of ([8|) yields 
~ (r(l/4) - r{l/A,{kxl)'^))/x']_ (F's are Gamma and 
incomplete Gamma functions). Last expression is 
plotted in fig. [1] function of kx']^ along with a 

constant line (see dashed line in the figure) of I /{GkE). 
The intersection point kx'[^ is the rescaled radius of the 
unique apparent horizon for the given values of E and k. 

It is evident that regardless from the value of GkE, 
which can be thought as a measure of the energy density 
on the transverse plane, a trapped horizon always exists: 
no mutter how dilute the energy is distributed on the 
transverse plane, a single BH will be formed. 



p can not drop faster than because the total energy of the 

collision is finite by assumption. 

In particular, this implies that (xp)' can not become zero in the 
interval (0, oo) because xp attains its unique (due to assumption 
(iii)-see below equation Js}) extremum at x = 0. 
Uniqueness can be shown using the fact of the previous footnote 
and working analogously to the corresponding proof of section 
1111 A 31 (see the proof in the relevant footnote of that section) . 



This case corresponds to distributions that near the 
origin behave as p{x±) ~ l/x^ -I- (sub-leading) This 
implies that the left hand side of ([5]) goes to a constant 
as a;'^ — > and goes to zero as x'^ — ?> oo and as a result, 
there is a single solution^^ x']_ to equation ([8]) but only 
for sufficiently large E. 

As an example, we consider the (everywhere posi- 
tive and integrable) distribution p ^ k'^e~^ ^^/{kx), 
which (for completeness) results to the shock (p ^ 
GEkxjL 2F2{{l/2,l/2};{3/2,3/2};-{kx:Lf). The left 
hand side of ^ yields ^Erf(fca;5^)^))/a;5_ where Erf is 
the error function. Previous expression is plotted in fig. 
[5] as a function of kx']_ along with lines of constant values 
of \/{GkE). The intersection point /cx'^, if it exists, is 
the rescaled radius of the unique apparent horizon. It is 
evident that a single BH will be formed if and only if the 
transverse energy density is large enough. 

An interesting limit of this example is the limit where 
— > oo. Then, the distribution becomes a delta function. 
This is precisely the case of Eardley and Giddings [l^ 
where they found that a BH is always formed. This limit 
is easily recovered by considering the ordering of limits 
limfe^oo[hm3;^^oErf(/cxj_)/a;_L], which is consistent with 
the ordering of limits that a delta function is defined 
and, which yields ~ limfc_>oo k. This implies that the in- 
tersection point of the left hand side of (|8]) with the ver- 
tical axis (see fig. [5]) is shifted to infinity. Equivalently, 
limfc^ooErf(fca;j^)/a::^ ^ l/x±^ which is the ,14] case and 



Since {xp)'\x=o = 0, it impHes that xp does not have any other 
extrema (due to assumption (iii)-see below equation 
Uniqueness can be shown using the fact of the previous footnote 
and working analogously to the corresponding proof of section 
IIII A 31 (sec the proof in the relevant footnote of that section) . 
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FIG. 3: The function (1 — )/y as a function oi y = kx±. 
This is an example where there exist two trapped horizons 
provided that the transverse energy distribution is dense 
enough. The dashed lines correspond to lines of constant 
E. The top and bottom lines imply no BH and two BHs for- 
mation respectively. The middle line suggests that there is a 
critical energy such that the large and the small BH merge. 



hence a BH is always formed in this limit. Effectively, 
this limiting case transforms into the previous class of 
geometries/ distributions we have already studied. 



3. Two black holes for sufficiently large energy 

Finally, the third class concerns distributions that be- 
have as p{x±^ ^ 1/a;" -l-(sub-leading) where n <1. Then 
the left hand side of ^ vanishes at a;^ = and since it 
also vanishes at = cx), it implies that given sufficiently 
large energy, there exist two roots for equation ([8]) result- 
ing to two BHs: a small and a large one. The possibility 
of more than two roots of ([8|) and hence more than two 
BHs is excluded as a consequence of assumption (iii) in 
section IHL^ ". 

As an example, we consider the (everywhere positive 
and integrable) distribution p ^ k'^e~^ The cor- 
responding shock is ^ EG (^Ei{k'^x\) — \og{k'^x\)) 
where the logarithm removes the logarithm of the expo- 
nential integral function Ei (when expanded at small ar- 
guments) and hence it removes the delta function source 
keeping only the gaussian one. 



Pf: The extrema of the left hand side of l(8j are the zeros of 
the function f{x) = px^ — pxdx and are located in the open 
interval (0, oo). We claim that / begins from zero (because of 
assumption (ii)), then for small arguments it increases at positive 
values (because at small x, p{x) increases slower than 1/x) and it 
eventually asymptotes to a negative value (because of assumption 
(ii) which in particular implies that p decays faster than 1/x-^ at 
large x.). This implies that / must have at least a (positive 
valued) maximum and at least one root. We claim this root is 
unique. To show this, it suffices to show that the maximum of / 
is the unique extremum, which means that /' must have a single 
root. Computing /' we obtain /' = x{xp)'. But /' ^ (xp)' 
has only one root in the open interval (0, oo) which is true by 
assumption. 




FIG. 4: The function x{xp)' as a function of x (in units of 
k) where p k^ (fcx cos'^(A:a;) exp[— O.lfca;]) . This distribution 
clearly violates condition (iii) and more than two apparent 
horizons are thus expected. 

The left hand side of dS]) yields (1 — e^'^^^i)/a;^ and is 
plotted in fig. [3] as a function of kx\^ along with lines of 
constant values of \/{GkE). The two intersection points 
/ca;'^, if they exist, are the two (rescaled) radius of the 
apparent horizons. Evidently, there are two trapped sur- 
faces, a small and a large one, if and only if there exist 
enough energy in the collision. In fact, the small trapped 
surface becomes smaller as E increases with opposite be- 
havior for the large one. What is more, is that there ex- 
ists a critical value of the energy where only one trapped 
horizon exists while for even smaller energies, trapped 
horizons cease to exist at all. These features point, by 
analogy, to the direction of an extremal charged BH and 
a naked singularity situation respectively (compare with 
[13] )■ We will return on this in section UlI Bl 



4- Many-many apparent horizons 

Relaxing the (iii) condition, which says that the ex- 
pression x{xp)' has a single root (other than maybe 
the root at a; = 0) and allowing more roots, gener- 
ally yields more trapped surfaces. As an example, we 
consider the positive and integrable distribution p ~ 
/c^ (fca; cos^(A:x) exp[— O.lfcx]) where k is the transverse 
scale associated with the width of p. According to fig. 21 
it violates (iii) condition. Hence, more BHs are expected. 
Indeed, fig. 2] shows the violation of condition (iii) while 
fig.'sOindicate explicitly the possibility of more than two 
apparent horizons may be created. 

An important observation is that at high enough ener- 
gies, the number of the resulting trapped horizons is (at 
most) two and, this independently on whether assump- 
tion (iii) applies. This statement holds because (the left 
hand side of ([5])) vanishes only at a; = cxd and depending 
on the behavior of p at the origin, occasionally at a; = 0. 
Lowering sufficiently the horizontal lines (see caption in 
fig. [5]) by increasing the energy, it will result in two inter- 
section points (at most) and thus into two trapped hori- 
zons exhibiting some universal behavior independently 
on the (intermediate) details of the p's. 
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FIG. 5: The left hand side of condition (|8]) for the distribu- 
tion p ~ fc^ (fca; cos^ (fca;) exp[— O.lfca;]) as a function of y = kx. 
The top figure shows the condition globally in the whole 
range. At large enough y's the left hand side of ((8| decays 
to zero. The bottom plot, focuses close to the origin. The 
number of BHs is found by the number of intersection points 
of the plot with horizontal lines, which correspond to some 
fixed energy. Clearly, the horizontal (dashed) line with value 
2.5 intersects the curve three times at the vicinity of y — 4.5 
and (at least) one at infinity yielding four possible trapped 
horizons. 

B. Thermodynamical implications 

We will discuss the small BHs of llll ASl with respect to 
their thermodynamical properties and thermodynamical 
stability. Most of the discussion of this part is based on 
assumptions and speculations and outlines several possi- 
ble scenarios. However, a resolution of the open possibil- 
ities would be achieved when the actual Einstein's equa- 
tions are solved and the problem is posed as an initial 
value problem in which case numerical relativity could 
say the final word. Our approach in this paper is simpler 
and inadequate in distinguishing the several ways of the 
evolution of the system, which will be mentioned below. 

We begin from the Hawking formula for the entropy of 
a BH 

(9) 

where Atrap is the area of the trapped surface, the factor 
of 2 (second equahty) comes from the two trapped sur- 
faces S+ and S- (see section (|ll|) while the dependence 
of the trapped horizon on the energy is explicitly given 
(third equality). In what follows we will assume S bh is 
proportional to Strap = S and drop any subscripts on 



Strap- Numerical works [1^ (in AdS backgrounds) indi- 
cate that the difference is an energy independent propor- 
tionality factor. 

We see two possible scenarios. 

1. The co-existence scenario where the small surface 
is inside the big one (co-centric) and they are both cre- 
ated at the same time. Then, the following possibilities 
exist: (a) Only the large BH mutters as the small one 
is hidden inside the large BH and its effects/presence 
will never be observed, (b) A related possibility is that 
the resulting black object will eventually evolve into a 
Reissner-Nordstrom (RN) BH. An RN BH has two co- 
eccentric horizons separating the BH in two regions. In 
this work, we find that the outer (inner) horizon increases 
(decreases) as the energy increases, similarly to the RN 
case. In addition there exists a critical value of the en- 
ergy (see middle dashed line in fig. [3]) where the two hori- 
zons coincide, resembling the extremal case in RN BH. 
The connection can be motivated further by defining the 
transverse scale of the width k to be Tq = = Q^G 
where Q a dimensionless parameter that can be dialed 
in independently on E. Hence, in our case as well, there 
also exist two length scales: rs = GE and = Q^G 
while the extremality occurs when rs/rq take a given 
minimum value. Decreasing the energy further, we find 
that a trapped surface is not formed. The corresponding 
picture in RN yields a naked singularity. 

2. The separately formed scenario: the two solutions 
are possible saddle points of the (classical) theory. How- 
ever since the large BH has a larger S, it has more prob- 
ability to be formed. But in principle either the one or 
the other will be formed but not both (simultaneously; 
classically at least). We are mentioning the following 
possibilities: (a) The small BH is unstable and will de- 
cay immediately, (b) This (the small) solution naively^''' 
yields negative temperatures because its size (and hence 
entropy) decreases as energy increases, hence, the small 
trapped horizon solution just doesn't make sense and it 
should be discarded. 



C. Universal behavior(s) 

In the high energy limit, the left hand side of ^ be- 
comes 1/xc because Xc becomes large and as a result, 

" xpdx — 1. Thus, Xc — >-~ GE with the dependence 
on k dropping out. This implies that in the high energy 
limit the following universal result^^, independent on the 



In the sense that at the instant of the collision, the system is far 
out of equilibrium and the discussion about temperature should 
be taken with caution. 

When two BHs exist, the universal behavior refers to the large 
BH. 
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FIG. 6: The (un-normalized) boosted Woods-Saxon distribu- 
tion plotted numerically as a function of the transverse coor- 
dinates (in arbitrary units). Evidently this distribution be- 
longs to the third class of distributions of this section and as 
a result, there expected two BHs if the energy is sufficiently 
large. 



FIG. 7: The (normalized) left hand side of ((Sj times a = 
0.5 fm plotted numerically as a function of the transverse 
coordinate (in units of fm). As expected from fig. ((6| and the 
discussion of subsection IIII A 31 two BHs may be formed given 
sufficiently large energy. The plot has a peak approximately 
at the point (13.5,0.024) and this point corresponds to the 
minimum energy where an apparent horizon could be created. 



transverse distributions, applies 

S G{Ekf ^k^s, EkG->l. (10) 

Another interesting point is to note that all the shocks 
in subsections IIII AlllTlIABI grow logarithmically at in- 
finity. In fact, this is valid for any (positive, integrable) 
distribution p. To see this, one must note that the deriva- 
tive of the shock is 4>' ~ l/a;_L /q ^ yp{y)dy and it be- 
comes (j)' ^ l/x± at large x±. This implies that in 1+3 
space-time dimensions, a (physical) distribution p, which 
creates a shock geometry such that it decays at infinity, 
does not exist. In 1 -I- n space time dimensions the shock 
decays as cj) ^ at large x±. As we will see, there 

also exists the AdS version behavior of the shock- wave at 
large arguments and, which specifies the stress tensor in 
the gauge theory side. 



D. An application for LHC 

An interesting distribution to investigate is the Woods- 
Saxon (phenomcnological) profile [l^, which seems to de- 
scribe the nuclei distributions quite successfully. Boost- 
ing to the infinite frame by sending the boost factor to 
infinity but keeping the total energy E fixed (at the order 
of the LHC scales), one obtains (see fig. H]) 



T++ = 



E 



dX3 



cV2 J-ao 1 + gVi^I7^r+<^I7^-R/>^ 



S{x+) (11) 



where a w 0.5 fm, R 7.4 fm is the Pb radius for the 
probes at LHC, c = -TrLi^^exp^R/ a)) fm^ = 13732 fm^ 
(Li3(x) is the poly logarithm function), which is chosen 
such that when (fTTj) is integrated on the whole space to 
give the total energy E. The \/2 is a consequence of our 
convention for the light-cone coordinates. 

In fig. ([7]) we plot the quantity on the left hand side 
of ([5]) by identifying k with 1 /a while both of the sides 



of the equation are dimensionless. As expected, it seems 
that two apparent horizons may exist, a small and a large 
one, given sufficiently large energy. On the other hand, 
the exact proportionality factor on the right hand side 
of ([8]) turns out to be 2\/2/(167r)a sa 0.03 fm and as a 
result, this side of the equation is approximately equal to 
0.03 fm/GE. Converting everything in Planck mass units 
(Afp —l/G) the right hand side of ^ becomes 



1.8 X W'^Mp/E. 



(12) 



Taking into account that, according to fig. [71 the trapped 
horizon with the least energy is formed when the left 
hand side of ([5]) is 0.024 and, equating this value with 
(|12p. it is deduced that the required energy is approx- 
imately 10^" orders of magnitude of the Planck mass, 
which is an extremely large energy. We conclude that, 
if this simplified model is not far from reality, no BHs 
will be seen at the LHC. However, there is a scenario dis- 
cussed in the conclusions, which may avoid the obstacle 
in creating BHs of such a large scale set by the Newton's 
constant in 4 dimensions. 

It is also interesting to note that any possible BHs 
found at the LHC (if the alternative scenario discussed 
in the conclusions is employed) could be of the RN type 
in the sense that they will consist of two co-eccentric 
horizons having a similar behavior with the energy as in 
the RN case. 



IV. AdS^ BACKGROUND BHS PRODUCTION 

Our investigation concerns head-on collisions of sym- 
metrical distributions that (may) yield 0(3) trapped ge- 
ometries. We classify the 0(3) symmetric bulk distri- 
butions in AdSs according to whether trapped surfaces 
may or may not be formed. This symmetrical case is cer- 
tainly not the most general one but it is convenient as 
it allows an analytical treatment of the problem and it 



8 



suffices in showing our claims. The presentation of this 
section, which (initially) follows [2l|, will be brief as the 
ideas are similar with those of previous section. We be- 
gin from the analogue of equations ([T]) and adapted 
to spherically symmetric shocks in AdS^ 

ds^ = 6(z)^( - 2dx+dx- +dx\+dz'^) + b{z)<j){q{x^ ,x'^ , z)) 



X S{x+)dx+'^, b{z) = 



L 



q = 



x\ + {z- zof 



4zzo 



(13) 



where L is the AdSs radius, z is the fifth dimension with 
the boundary CFT residing at z = 0, q is the chordal 
distance in AdSa, which makes the 0(3) symmetry man- 
ifest while the role of zq is analyzed below. This ansatz, 
as before, satisfies trivially all the Einstein's equations 
except the (H — h) component, which yields 



i?4 



=> -l/2S{x+) (D^L - (f> = 87rG(z/L) J++ 

= 87:EGp{q)d{x+), J++ = E{L/z)p{q)5{x+) (14) 
where the □.l operator is with respect to the AdS'i metric 
ds^ = £2 [dq^li^qi^q + 1)) + 4^(1 + q)d^-^^ , (15) 

G is the Newton's constant, which in this section con- 
cerns five dimensions, E is proportional to the total bulk 
energy of the distribution, zq is the location of the cen- 
ter of the bulk distribution and, which according to the 
AdS/CFT dictionary (see [HI, HI] for instance), is pro- 
portional to the width of the gauge theory stress-energy 
tensor (see also below). Both of the quantities (j) and p are 
functions of q and they have dimensions of length to the 
first and the minus third power respectively. Dropping 
any positive dimensionless constants but keeping track 
of relative signs, previous equation yields the following 
differential equation for (p 

q{l + q)<j)" + 3/2(1 + 2q)<j)' - 30 - -EGL^p 

= -EG/L^zqp (16) 

where p = L'^/zoP as defined is dimensionless while cj) 
has dimensions of length as it should. The choice of 
the overall coefficient for the bulk source comes about 
as follows: we need the field (j> to expand at large q as 
1/q'^ in order to ensure consistency with the AdS/CFT 
dictionary. Then, taking into account (|16p and the fact 
that T++ ~ (L^/G) X l/{zq)^\^=oSix+) where T++ is 
the gauge theory stress-tensor and assuming that we 
fix the integral of on the transverse plane to yield 
E^^, we conclude that the overall dimension- full coef- 



Allowing us to identify E with the total collision energy in the 
gauge theory side. According to II14I I. E is (at the same time) 
proportional to the total bulk energy of the collision from the 
5-dimensional point of view. 



ficient of </> is EG/LP'zq. Indeed, such a choice results 
^ E , ■iy^ Six'^) which integrates to E. 



The solution to equation (|T6|) can take the integral 
form 



E{G/L')zo{{l + 2q)[Gi+ r ai{y) r a2{x)dxdy 

Jqo Jo 



C.|^^il±Mi±^+8(l + 2,) 



1 



"^^'^^ (l + 2g)2(q(l + g))3/2' 
a2{q) = {l + 2q)^q{l + q)p{q) 



(17) 



where Ci(go) and C2 — J a2{yo)dyo\y^^Q are constants 
chosen such that to remove the (l-|-2q) and l/^/q terms 
of (j> coming from the two lower limits of the integrations. 
This choice of Ci yields (as mentioned above) </) ~ 1 /q'^ 
at large q and ensures that the geometry, according to the 
AdS dictionary, corresponds to a boundary CFT theory. 
C2 on the other hand, should be chosen such that to can- 
cel the l/-y/g behavior of the integral at small q because 
such a behavior induces an additional (to p) ^(g)-source^^ 
on the right hand side of ([T6l) (see [2l|). In particular, 
for distributions satisfying a2{yo)dyo — >■ as y — ^ 0, 
the coefficient C2 = 0. 

A few remarks are in order: (a) the first integral in 
()16p involving 02 not only it should be finite as the up- 
per limit tends to infinity but it should yield a positive 
constant. The reason is because, as we will see soon 
(see equation ()19|) ). it represents the bulk energy inside 
a ball of radius y sitting in AdSa and this energy should 
be finite. This obviously occurs for a subset of choices 
of positive p's; we call these physical p's. (b) In par- 
ticular, remark (a) implies that in the large upper limit 
expansion, this integral has a zeroth order term while 
the rest terms are sub-leading (involving negative pow- 
ers), (c) The second integral involving oi when convo- 
luted with negative powers and then expanded in the 
large upper (q) limit, yields powers that drop faster than 
l/q'^- Hence, these terms do not contribute to the gauge 
theory stress-tensor. Convoluting with a zeroth power on 
the other hand, yields an 1 /q'^ fall-off and this term does 
contribute to the CFT stress-tensor, (d) Previous three 
remarks imply that physical sources yield a proper map- 



Pf: In order 

grate (DAdSa 



with / 



to see this one should, using l|15p . inte- 
— 3/L^)f inside a small sphere of radius qo 
and finally take the limit qo — > 0: 




0{^) - 



-1/2 /s-?!! 

where g™'' the induced metric of 33 _ y^ia^g 
S . q = qo, a unit normal vector on the same surface while 
after first equality, Gauss law has been employed. 
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ping of the geometrical description to the gauge theory 
description^* . 



A. General classification 



The consequences of these three assumptions are simi- 
lar as in flat space (see discussion below In particu- 
lar, condition (iii) ensures that no more than two trapped 
surfaces may be created^". 



The analogue solution of dH) , to equation ^ with the 
condition ipiQc) = satisfied and where the curve C : q — 
q'^ defines the trapped horizon, is 



V'('7) = (l>{q) - 



4>{qc) 

l + 2q. 



(l + 2<z). 



(18) 



The derivative condition on the trapped surface in this 
case is g'^'^ip ~ 1 where g'^'^ = q{l + q)/L is the (qq) com- 



ponent of the inverse metric of AdS^. 
p7| and (Ull) yields 



Then, employing 



il + 2qc)qc{l + qc) ' 



ZQ GE 
k 

'GE' ^' 



1/zo (19) 



where ~ l/^O) as mentioned, is a transverse scale asso- 
ciated with the width of the gauge theory stress-tensor 
in the CFT. Equation ((T^ is the same as the condition 
derived in [2l'| and has an analogous form and physical 
meaning as condition ([5]) with one crucial difference: the 
right hand sides of the two conditions differ by the fact 
that in there exists the combination Exk while in 
there exists the ratio^^ k/E. This is a very important fact 
with potential consequences in heavy ion collisions and, 
which we will return in section flV CI 

The comments below condition (jS]) apply also here with 
minor modifications and so we will be brief in the subse- 
quent presentation. We have therefore, all the machinery 
and the ideas in classifying the 0(3) invariant distribu- 
tions in AdS according to whether they produce one or 
more trapped surfaces during a collision. The distribu- 
tion is again assumed (i) everywhere positive, (ii) inte- 
grable on the AdSs subspace and (iii) that the quantity 
yjq{l + 2q){l + 2q)p{q)/{l + Qq + &q^) has at most one 



1. Always a single black hole for any energy 

This case corresponds to distributions that at small q 
behave as p ^ l/g"-|-(sub-leading) where 3/2 > n > 1/2. 
The upper bound comes by demanding finite total energy 
while the lower bound comes by demanding the left hand 
side of (HI) to diverge as qc 0. 

As an example of such a distribution one may 
consider p ^ 1/['7(1 + 2(7)(1 -I- qY], which as can 
be checked, creates the shock (f) ^ —EGzq/ L'^{1 + 
2q) (1/(1 + g) + 4/(1 + 2q) + log[g(l + 29)7(16(1 + qf)])- 
At large q this shock-geometry behaves as l/?"^, which 
implies the stress-tensor in gauge theory behaves as 
E/ {x\_ + ZqY and is also positive as should. It is pointed 
that the relative signs work out^^ and that zq is the trans- 
verse width of the gauge theory stress-tensor. The left 
hand side of (dH) behaves as 1/ {y/q{l + qf/'^{l + 2q)) 
and the corresponding plot topologically looks like 
fig. [1] showing that a single trapped surface is always 
formed for any energy. Evidently, at large E one obtains 
ql - E/k. 



2. Marginal case: one black hole for sufficiently large 
energy 

This case corresponds to distributions that at 
small q behave as p ^ 1/^-1- (sub-leading). As 
an example of such a distribution one may con- 
sider p - 1/[(1 + 2qf^{l ^ 

be checked, creates the shock 



4V9(1 

-(1 + 



+ g)(l + 2q){t&n-^ 
8<z(l + g))log((l + q)/(l 




g)'^/^], which as can 
EGza/L^ (2q - 
log(2) - 7r/4) 

2g)) )/V9(l + 9)- At 



maximum. 



The l/q^ decay of <f> at large q can be justified rigorously as 
follows. The term in II17I I which is proportional to C2 decays 
as and so we ignore it from the further argument. Divid- 

ing the remaining terms with (l-|-2q), taking the derivative in 
both sides of the resulting equation and taking the large q limit, 
taking into account that the integral involving 02 is finite by as- 
sumption (see I I19II and assumption (ii) in section [IV At . yields 
{4>/{l + 2ij))' ~ l/q^ where <j > 1. Integrating (0/(1 + 2q))' (at 
large q) yields ~ l/g^ up to 1 + 29 terms coming from the in- 
tegration constants and which are cancelled by a suitable choice 
of Ci in flTll . 

The dimensionless parameter G^/L^ is removed from the discus- 
sion as, in a bottom-up approach, it will be chosen such that it 
fits data [U] . 



Working analogously as in the fiat space, one finds that a con- 
dition for having two trapped surfaces at most, reduces in de- 
manding that the function g{q) = f /b where / = \/q{l + 2q)(l + 
2q)p(q) and b = (q{l + q){l + 2q))' = {l + Qq + Qq^) has (at most) 
a single maximum at some point q = qo. Remark: The condition 
that the quantity / = ^ q{l + 2g)(l -|- 2q)p has at most a single 
maximum is more physically motivated by analogy with the fiat 
space and in the view of the integrand of l|19|l . It is possible 
to make the two conditions applicable simultaneously by impos- 
ing extra conditions on /. For instance one could demand that 
(log(gfe))' = (log(/))' has only one root (at most) in the interval 
{qo,oo) where qo is the (only possible) maximum of g. In fact, 
the two conditions happened to coincide for all the cases that we 
have checked. 

Such that both, the (dual) gauge theory stress-tensor and the 
bulk tensor are positive. 
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large q the shock decays as l/(?^, which is in accordance 
with the AdS dictionary while, the left hand side of 
(HI becomes ~ log[(l + 2q)/{l + q)]/{q{l + q)il + 2q)), 
which vanishes as g — >■ oo and it is finite when q = 0. 
Topologically, this plot is like that one of fig. [2] showing 
that only a single trapped surface may be formed for 
sufficiently large energy. In addition, at large E one 
again obtains q^ ^ E/k as in IIV A ll 



3. Two black holes for sujjiaently large energy 

This case corresponds to distributions that at small 
q behave as p ~ 1/(7"+ (sub- leading) where n < 1/2. 
As an example of such a distribution one may consider 
the n = —1/2 and in particular the distribution p ^ 
^/ ((1 + g)5/2(l + 2(j)2), which as can be checked, cre- 
ates the shock (j> EGzo/L^(3 + 2q{30 + q{61 + 32q))- 

36(1 + q)^q{l + q)il + 2q) tan"! (^^/(l + ?)) + 6(1 + 

q)il + 8q{l + q)) tanh-^ (q/(2 + 3g)) )/ (3^9(1 + q?/^) + 

EGzo/L^ ((l+8g+8q2)/ y^(TTg)) +EGzo/L^ ((4/3(5-1- 

log(8)) — 37r) (1 -I- 2q)j . The second and first terms from 

the end subtract the delta function source at small q and 
ensure a 1/q^ (in order to have a proper CFT mapping) 
fall-off at large q respectively. Remarkably, this shock- 
geometry is everywhere finite, has finite first order deriva- 
tives but irregular second and higher order derivatives at 
q = 0. 

The left hand side of ([19]) becomes ^ ('^ 

(1 + g) log [(1 + q)/{l + 2q)] ) / {q{l + q)^l + 2(7)), which 

vanishes as g and as g ^ cx). Topologically, this plot 
is like that one of fig. [3] showing the possibility of cre- 
ation not of just one but of two trapped surfaces, a small 
and a large one, for sufficiently large energy. Finally, the 
behavior of qc with E in the high energy limit is the same 
as in previous cases, namely as q^ ~ E/k. 

One should bare in mind the a naly tical work of [s^l 
and the numerical results of [HI, l52j | concerning AdSs 
background. In these works, a BH eventually, after mul- 
tiple re-scatterings with the boundary of the AdS, will 
be formed. However, these works have a different set-up 
from the current one. They use a global covering of the 
AdSs space, the corresponding bulk fields are self inter- 
acting in the sense that they obey their own equation 
of motion and most importantly, these are exact confor- 
mally invariant set-ups. In the contrary, in our case the 
presence of the energy E and the transverse width k of 
the distribution breaks conformality (softly). Hence the 
fact that a BH is not always formed is a possible sce- 
nario as supported by our results and also from other 
works [53| . 



B. Thermodynamical implications 

For AdS backgrounds we know that the AdS- 
Schwarzschild large BH is thermodynamically stable and 
the small one is unstable. Thus, it seems that all the 
trapped horizons that we find and, which eventually will 
form BHs, will be probably thermodynamically stable ex- 
cept maybe from the small horizons of subsection IIV A 3l 
(in a non-coexisting scenario see section UlI A3|) . In the 
coexisting scenario, where both, the large trapped sur- 
face surrounds the small one, it is possible that the sys- 
tem will evolve into an AdS-RN BH where, as in the fiat 
space, the third required scale (in addition to the AdSs 
radius and the Schwarzschild scale due to E) is probably 
the transverse width of the distribution. 



C. Implications for heavy ion collisions and QGP 

At RHIC and now at LHC the QGP production is 
expected that is formed during heavy ion collisions. 
AdS/CFT has been proven a powerful tool in modeling 
these kind of collisions (H-HS El, El UMl and predict- 
ing the QGP formation, which in the AdS/CFT context, 
implies a BH creation. However, the existing literature 
has missed one important feature of QGP observed ex- 
perimentally and, which is according to our intuition: 
QGP, even for very central collisions, is formed for suffi- 
ciently large energies, much larger than Aqcd- However, 
the existing works have found that the criterion in form- 
ing a BH shares qualitatively the features of the plot in 
fig. [TJ a BH is formed for any small energy (as long as 
the collision is central). 

Identifying k with Aqcd in (I19p is now tempting. Ac- 
cording to (the analogues of) fig.'s[T] and [5] (adapted for 
the AdSs background case), it is deduced that at high en- 
ergies such that E ^ k = Aqc d , there are more chances 
to get a BH and hence to produce QGP. Opposite con- 
clusions apply in the <C fc limit. This is exactly what it 
is expected/observed in realistic heavy ion collisions and 
it is pleasing that it can be modeled geometrically. 

In this work, we have shown that one can built geome- 
tries that both, map properly onto a CFT and, which 
when applied in heavy ion collisions, they capture the 
prescribed feature of QGP formation. We assumed an 
0(3) symmetry for simplicity, which allowed us to de- 
rive our conclusions using purely analytical methods. 
This symmetry, which is not necessarily realistic (in real 
heavy ion collisions; see [76j however) together with (ap- 
proximate) conformal invariance, force the gauge theory 
stress-tensor of the colliding matter to behave as a power 
law and in particular as ~ Ezq/{x']_ + z^)^S{x~^). In 
addition, as we have seen from all the examples we have 
studied, in the high energy limit qc - {E/kf/'^. This im- 
plies that total multiplicities N grow as S ^ {s/k'^Y^^ 
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where s is the cm. energy of the heavy ions beam^^ . This 
is a general feature of the 0(3) symmetry and yields a 
universal behavior^-^ for S{E) similar to the flat space 
of section IIII CI Such a universal high energy limit be- 
havior predicts a slightly larger growth of multiplicities 
compared to the data {N ^ s^/'*)^'*. In order to see this 
universal result, one has to take the high energy limit and 
hence the large Qc limit of taking into account that 
the distribution p is integrable (by physical assumption). 
Hence the numerator of the left hand side of this equation 
approaches a constant and the equation yields ~ E/k. 
The final step is to compute S from the adaption of Q 
for AdSs with the subspace measure (that of AdSa). At 



large qc it yields S ~ qc^^ and hence S ^ N ^ [s/k'^) 



2U/3 



Attempting to fit the RHIC and the LHC data by incorpo- 
rating weak coupling physics: a phenomenological approach 

It is possible that the scale k could depend on the en- 
ergy. Instead of identifying the transverse scale of the 
gauge theory stress-tensor k with Aqcd one could iden- 
tify it with Qs{E) that is with the saturation scale (78j . 
Such an identification is physically motivated as Qs is the 
typical transverse momentum of the gluon field (which 
dominates over the fermions) in a highly boosted nucleus. 
Furthermore, taking k — Qs, incorporates, in a sense, the 
weak coupling physics that also participate in heavy ions 
collisions. Taking into account that perturbative calcula- 
tions [li-lsa yield Qs ~ Aqcd(v^)^, A e [0.1,0.19] with 
A = 0.19 [S^l referring to nucleus-nucleus collisions and 
s is measured in GeV, it is deduced that^^ 

N^{^]'^'^ ^(./A2 )0.27_ (20) 



Such a dependence of on s, which is very close to 
the commonly accepted power N ^ s^/* (see [l^ and 
references therein), may fit the Phobos [sJl and the first 
ALICE [s^ data satisfactory within the error bars, up 
to an overall energy independent proportionality factor 
(see fig. [8] and 23] for a different but related approach). 
This result is universal, independent on the details of the 
transverse distributions, as long as the high energy limit 
is assumed. This implies that E — 1/2-^5 k — Qs 



The fact that Af ~ S is based on phenomenological assumptions 
and it is not straightforward. The fact that S ^ s^^^ will be 
explained shortly. We consult the reader to see [T9l423ll . 
Departing from the 0(3) symmetry, different (than the power 
law) gauge theory stress-tensors become possible and one may 
find a model that is closer to more realistic nuclear profiles. To 
this end, multiplicities are generally going to be different, 
bee [23|, [77] for a more involved treatment where a fitting with 
RHIC and LHC data as well as future predictions are achieved. 
Equation II20I I may have significant corrections due to both, 
higher order corrections and also running coupling effects of A 
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FIG. 8: The PHOBOS [8J[ and the preliminary ALICE [85t[ 
(right point) data including the error bars. The vertical axis 
denotes the total multiplicity while the horizontal denotes the 
nuclon-nucleon cm. energy meassured in GeV. The far right 
(ALICE) point has been adjusted from [s^ by a factor of 
2Apb/375, where Apb is the atomic number of Pb, in order 
to compensate the non-exact centrality aspect of the process. 
The dashed line is our theoretical prediction, (formula (|20p ') 
iV = 295{^/GeVf'^^ where the number 295 is the only 
(energy independent) fitting factor. 



should apply, which means (y/sf-AqcD)^ S> 1 should 
hold and, which is a consistent inequality for the energies 
at RHIC and at LHC. 



Mapping the gravity evolution to the evolution of the 
produced medium m heavy ions 

A possible evolution of the system in gravity and its im- 
plications in heavy ions motivated from the gravity anal- 
ysis of section IIV A3I is the following. After the shocks 
collide, two co-eccentric 3-dimensional black disks, the 
trapped surface, is formed. This surface has (almost) 
a zero depth in the longitudinal direction. The projec- 
tion of this configuration in i?^''^ is 2-dimensional disk, 
that is a 2-dimensional plasma. This applies immedi- 
ately after the collision. Then, the trapped horizon will 
isotropize due to the huge longitudinal pressure gradient 
in the center of the black object and it will evolve into a 
(RN) BH; exactly as it is expected in heavy ions (86l - [88| 
where the resulting thermal medium becomes isotropic. 
Such an RN-like configuration could capture the pres- 
ence of a chemical potential (2^] in QGP. Since the in- 
ner horizon of an RN BH has different thermodynami- 
cal properties than the intermediate shell extending to 
the outer horizon [89|], this system could model a hotter 
central region and a thermodynamically unstable periph- 
ery |3ll - [34| during the medium's last stages of evolution. 
Finally, the BH will undergo through a Hawking evapo- 
ration, which implies a hadronic phase transition in the 
heavy ions picture. This difference in the thermodynam- 
ical behavior of the inner-outer region is qualitatively in 
accordance with the hydrodynamical simulations [90l49^ 
of the QGP evolution. 
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V. SUMMARY, CONCLUSSIONS AND 
FUTURE DIRECTIONS 

In this work we have studied the problem of trapped 
surface(s) formation resuhed by colhding extended 
"spherically" symmetric, physical sources on the trans- 
verse plane and boosted to the speed of light. The total 
energy of the collision has been assumed finite and the 
impact parameter zero. The investigation has taken place 
in 4 dimensional flat and AdSs backgrounds. 

Our conclusions are summarized as follows. 

I. For both backgrounds the spherically symmetric dis- 
tributions have been classified according to: (i) always a 
single apparent horizon for any E (ii) a single apparent 
horizon for sufficiently large E (iii) two apparent horizons 
for sufficiently large E. Last category suggests that small 
and large BHs may be formed even in flat backgrounds, 
which, at the best of our knowledge, is a new result. 

II. We found that the creation of small BHs in both 
types of backgrounds is possible. In section UlIB I we have 
presented the possible scenarios for the flat backgrounds 
with similar possibilities for the AdS backgrounds. We 
find particularly attractive the suggestion that the result- 
ing configuration will finally equilibrate into an RN BH. 
The role of the scale multiplying the charge in the RN 
BH should be played by the transverse scale of the distri- 
bution. In particular, considering AdSs backgrounds, the 
RN scenario would explain the appearance of a chemical 
potential during heavy ion collisions. Such a configu- 
ration could model a hotter core (using the inner hori- 
zon) having different thermodynamical properties than 
the outer region (the shell between the inner and the 
outer horizon) as it happens in realistic collisions, accord- 
ing to hydrodynamical simulations [90l - l9^ . A detailed 
scenario based on the RN BH explanation, with empha- 
sis in the implications of the QGP evolution, is proposed 
at the end of section lTV CI In any case, all scenarios worth 
further investigation while possibly only numerical rela- 
tivity may answer these sort of questions reliably. 

III. There is a universal behavior independent on the 
colliding bulk distributions, which have a transverse scale 
k that generally depends on E. In the high energy limit, 
one obtains S ~ (kE)^ and S ~ {E/kf/^ for flat and 
AdSs backgrounds respectively for any collided distri- 
bution obeying standard assumptions found in sections 
mil IIVI This applies even when assumption (iii) (see sec- 
tions mil and lIVp is relaxed. In the introduction, it was 
claimed that this result was independent on the presence 
of any extra dimensions and on the impact parameter b 
as long as the energy is sufficiently high. For the extra 
dimensions, this was shown in [93j in higher dimensions. 
The idea is clear: at very large energies, the resulting 
black object occupies the whole compact directions and 
becomes independent on them (see for instance (j23|) for 



large xx),. An analogous idea motivated by the exact re- 
sult of [14| , applies for having a non-zero b which results 
in a non-zero angular momentum and yields an ellipti- 
cal trapped horizon. However, as Eb — > cx), the surface 
becomes (more) circular. 

IV. A phenomenological model for describing possible 
(4-dimensional) BH creation at the LHC has been con- 
structed using the Woods-saxon profile as the colliding 
matter. This profile results to two trapped horizons given 
sufficiently large energy (see fig. [7]). Substituting the 
numbers, it is found that the minimum energy required 
is of the order of 10^° of Mp, pointing to the direction 
that no BHs will ever be seen at the LHC. This result is 
consistent with the fact that the SM interactions inside 
the world as we know it today, dominate over the gravi- 
tational ones for energies way below Mp. However, as it 
is discussed below, there might be a way around. 

V. Classes (ii) and (iii) of the first conclusion (see 
above) show the possibility in building geometries that 
both, mimic heavy ion collisions and, which predict QGP 
formation only for sufficiently large energies compared to 
the transverse scale involved and, which is tempted to be 
interpreted as the Aqcd scale. This fact applies (even) 
for perfectly central collisions. In addition, the 0(3) ge- 
ometries seem to predict that the total multiplicities at 
high energies behave as (s/k^Y^^ independently on the 
details of the profiles of the bulk distributions. 



VI. In the AdSs space, the transverse scale of the bulk 
distribution coincides with the transverse width of the 
gauge theory stress tensor before the collision and, which 
in turn, it can be thought as a nucleus moving with the 
speed of li ght. In the color glass condensate language 
(CGC) [zl [ill tliis width is the saturation scale Qs- It 
is of the order of GeV and it depends on the energy in 
a fashion which is computed perturbatively (95j . Identi- 
fying k with Qs is not only a plausible assumption but 
it also incorporates in a sense, the weak coupling physics 
in the large coupling treatment that AdS/CFT relies on. 
Such an identification yields the result that total mul- 
tiplicities grow as iV ~ (s/Aqcp)°'^^ which seems to 
describe satisfactory the PHOBOS [84] and the ALICE 
(ssl data at RHIC and at LHC respectively. The corre- 
sponding plot is depicted in figure [8] 



VII. A number of investigations have been postponed 
for future work. The non-zero impact parameter case cor- 
responds to more realistic situations, especially in heavy 
ion colhsions (i.e. elliptic flow [13]) and should be con- 
sidered. However, we expect that it will not change the 
results in the high energy limit (see conclusion IV and 
the results of [ll |4l| in the limit G£'/6 > 1). 
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Introducing extra dimensions and BHs at the Te V scale 



VIII. A more interesting new direction is the general- 
ization to more dimensions. Although a straightforward 
task, such a generalization will turn out useful in tak- 
ing another step further: the inclusion of extra, flat di- 
mensions for flat backgrounds, which is a valid scenario 
for the BHs production in the accelerators as it could 
lower the Planck mass and make the production of BHs 
feasible even at TeV scales. Related to this, we make 
the following hypothesis followed by a conjecture based 
purely on dimensional analysis. We require the creation 
of BHs at the LHC for collision energies of a few TeV's in 
the presence of extra dimensions having magnitude i?^, 
i = 1,2, ...n. We also assume that the colliding matter 
from a 4 dimensional point of view, belongs to the class 
of distributions of subsection IIII A 3l which are the more 
realistic ones and, which allow a BH creation only if the 
energy is above a threshold that typically is of the order 
of the Planck mass (Afp) or beyond. We then conjecture 
that creating BHs at lower scales than Mp is possible if 
the width = l/k^- of the matter in the extra directions 
is much smaller than the size Ri of the compact dimcn- 
sion(s). In other words BHs could be created at the LHC 
when the matter is sufficiently localized in the extra di- 
mensions such that kn.Ri 3> 1 is satisfied. Argument: 
since the Newton's constant when n compact dimensions 
are present, is G'4+„ — GiRi...Rm then the dimension- 
less coupling controlling the creation of the BHs should 
be 



1 



1 



M, 



,t2 



1 



Gi+n EkkR^...kR^ Ek {kR^Ri)...{kR^Rn) 

G4 ~ Mp ^ 

(21) 



'ff- {kn,Ri)...{kRM 



Ek 



where E and k are the total energy and the trans- 
verse width from the 4 dimensional point of view while 
Mf.ff is an effective Planck mass. Hence, for w ^ 1 
and taking into account that Mp/{Ek) ^ 1, then the 
R/fki = kR.Ri ^ 1 condition follows. This condition 
should had been expected. The reason being when this 
strong inequality is satisfied, the matter is spread along 
the 3-brane but it is localized (almost like a delta func- 
tion) in the extra directions. But a delta function distri- 
bution belongs to the first class (see subsection IIII A l( ) 
of distributions, which yield a BH for any small collision 
energy. Therefore, it seems that the localization of the 
matter in the extra dimensions, compensates the (realis- 
tic) spread in the two (transverse) extended dimensions 
allowing the possibility in creating BHs at lower scales. 

Taking E ^ IQ TeV and fc - 1 GeV equation ^ 
takes the suggestive form 



^ = TTT-^ ^Iw = Ek : 



where Mew is of the order of the electroweak scale. In 
particular, for n = 2, Ri = R2 = R ^ 0.1 mm as in 
[96l [qtI ] and demanding w '-^ 1 it is deduced that kn^ = 
kR^ ^ kR lO^Vrnm = lO^o/fm- 10^° GeV and so 
rfe„ = l/kR ^ 10-30 fm. 

It is remarked that such a (small) spread (osMlOl in 
the compact directions not only is allowed, as it does not 
violate anything fundamental, but is also essential from 
a quantum uncertainty point of view and thinking about 
the problem semi-classically. It is experimentally allowed 
for SM interactions to penetrate the extra dimensions as 
long as the corresponding penetration length^^ is much 
smaller than the resolution of the accelerators .The en- 
ergy required to resolve such small widths are of the order 
of the Planck mass or even (as in the example above) way 
beyond. Hence such a small penetration depth is practi- 
cally unobservable. 

As a candidate distribution, which could model the 
problem of not creating BHs in the absence of extra di- 
mensions at a given energy but with the possibility of 
BHs creation in the presence of extra dimensions at that 
same energy, one could use the distribution 



fc^/c^ exp [—k-^ 



(sin2(0i/2)+sin2(02/2))] 



47r3 (4i?2)e-'=H^'/2(ifc|i?2) 



(23) 



This distribution, which is normalized to unity and where 
Iq is a modified Bessel function, has the following fea- 
tures: (i) It is periodic with respect to the angles 9i and 
02 that parametrize the extra dimensions and it peaks in 
the vicinity of 9i = O2 = 2mT, n G Z where SM resides^'', 
(ii) Ignoring the extra dimensions, it belongs to the class 
of distributions of llll A3l where two BHs will be created 
for sufficiently large energy. This can be seen if one inte- 
grates out the compact coordinates. Then, the resulting 
curve is a gaussian and topologically it is what one ex- 
pects for realistic distributions in the 4 dimensional point 
of view (it looks like the one of figH]). This is what hap- 
pens at large energies. The trapped surface becomes very 
large and covers the whole compact directions [11] (it 
becomes a black string) and the distribution effectively 
becomes 2 dimensional depending only on the transverse 
coordinate x±. (iii) When the collision energy is suffi- 
ciently small, the trapped surface is too small to resolve 
the finiteness of the extra dimensions (see |93i] for a de- 
tailed analysis on this). Thus, when the the distribution 
is expanded (at small angles and small x±), yields a ro- 
tationally invariant expression if in addition, kji — 4fc is 
assumed. Then, the problem in this energy regime can be 



(100 GeF)^ 



(22) 



From a 4 dimensional point view such a characteristic length can 
be taken as the proton's radius. In the extra dimensions scenario, 
the length is the width of the distribution. 

Modulo a small but unobservable penetration inside the compact 
directions. 
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solved exactly and this would provide cross-checks^* of 
a more thorough (numerical) analysis, (iv) Plotting the 
distribution for fixed x±, it looks very localized; almost 
like a delta function. In fact, in the limit k^R — cxi it 
becomes a delta function while the denominator of ((23)) 
becomes kfiR independent. 

This distribution, from a 4 dimensional point of view 
(so ignoring the angles), creates BHs when Mp^ /{Ek) ^ 
1 is fulfilled. The question is whether the introduc- 
tion of the extra dimensions could yield BHs even when 
Mp^/{Ek) ^ 1 applies. The problem of solving the 
underlying trapped surface equations could be attacked 
numerically and answer quantitatively these questions. 
However, the intuition we obtained suggests that it is 
possible to have BHs in the Mp^ /{Ek) ^ 1 regime as the 
presence of the extra dimensions could lower Mp down 
to Meff ^ Mew yielding a BH at accessible collision 
energies. Hence, under the assumptions that extra di- 
mensions do exist and that the matter is dense enough 
in those directions, implies that LHC is maybe not that 
safe after all. 

Along with the high energy limit-see (ii) above. 
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